unc  FILE  COPY  AO  All  8710 


& 


urna  o»-  ha  via.  ta^uasN 


Contract  HoOOM-MJ-a-OOl  7 


Tui  Hu.  Mi  :•>•»- 1  la 


tCCMUCAt  AEPOMT  MO.  j  4 


SHuliUon  of  C4t tm  KUcu  in  l.l«ctiuai.4lytK4l  UiwiiH-nu  by 
Oitiu^oiMl  Collocation.  fill  III.  Application  iu  anuru^tu- 
Mtru  U^riWnU. 


Ihrnid  ^uci  «nd  Stanley  fun* 


Prepared  tor  hAhcitian  In 


KieCtrOChimica  Ac t A 


UtPOKT  DOCUMENTATION  PAGE 


,r  c  Kl  «U  lUWIUrilONl 

_ |iCH«*a  rn-tPi  tnw.  t  f> mi 

M^l  MClUttK  HU  l  M(>m  «  I  J  («UlM  m«<i  n 


Ik.  im  o>  Mi*o*t  e  ritMO  Ckvi  m  « 
Technical  Report#  l4 


I  *  ei»»  OMMinb  non.  Mr«ai 


I ».  ubi**ci  om  c*ui 


N00014-  82-0-0017 


_ L! _ WP'4  /Li 

»  Emulation  of  Edge  Clfit.ll  in 

Elwctroanulytii al  Experiments .  Hart  111. 
Application  to  ChronoaMpeiometc ic  Expci  Manti. 


fesrnd  Speiser  and  Stanley  Hons 


•  M  MOWmCOKimIMN*  NUC  Ml  iMMU  «MiVwM  VmV'wiSmMM  ‘  '4I* 

Department  of  Chemistry  task  No.  NR  359-718 

University  of  Alberta 

Frteir.nl  on,  Alt.erta.  Canada  Its  2G2 _ 

•  i.  (ONiMuwbanictawt  uomxmoi  u.  Meoar  kill 

Office  of  Naval  Research  _ _ 

Chemistry  Program  -  Chemistry  Code  472  a  "•“y,4!0'  '**“ 

_ Arl i naton.  Vir.nnia  ???1? _ _ *t= - 

I*  HomilOhi.C  scotl  cut  •  «•—  II  HCtOlff  Ulll  faMewMl 

Unclassified 

tt*"eic7  S  i  w  ic  *  ROSS SSKaiSS" 

ItalMtl 

M  fair«HiliM  UiUauT  H m> *«»— «j 

This  dociment  has  been  approved  for  public  release  and  sale;  Its  distribution 
w»l  ini  ted. 


August  1(1.  !?<U 


Uni  varsity  of  Alberta 
Department  of  Chemistry 
Montoo  .Albarta .  Canada 
T6G  2C2 


August  16,19*2 


((•production  in  who  la  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


This  document  has  been  approved  for 
And  salat  its  distribution  is  i 


or  public  rels«^i^^\  f  . 
s  unlimit 

W  ic^e 

W  VA 


AitexTw  HiUNUi  fa  awn  ■ 


I  >•  hlWUMUNIMIU 


Electrochemistry,  Chronoemperametry,  Simulation,  Orthogonal  collocation. 


ee«ie*CI  |Cmhm«wmm(W>  WMSMMy  ■eiawi»«rua«ep»,» 

^Comparisons  ot  simulated  results  with  theory  tor  chronoOMperometric  systems  is 
7  presented.  Aesults  foi  complicated  systems  are  obtained  with  high  accuracy  in 
very  short  periods  of  time.  ^ 


00  /(L,1*)  U7J  10.1.0*  S»  IMVttUS 

t/M  MSI 


Unf  1  .iislfifil 

MCWNM  U a.asnwm 


SIMULATION  OF  EDGE  EPFECTS  IN  ELECTROANALYTICAL  EXPERIMENTS 
BY  ORTHOGONAL  COLLOCATION.  PART  }  •  APPLICATION  TO 
CHRONOANPERONFTRIC  EXPERIMENTS 


hrnd  Speiser  and  Stanley  Pons 


Department  of  Chemistry 
University  of  Alberta 
Edmonton ,  Alberta,  Canada 
TAG  2G2 


Part  1;  see  ref.  Ill 


ABSTRACT 


Hie  theory  of  two -dimens  Iona 1  collocation  has  been  employed 
to  simulate  a  model  describing  chronoamperometr ic  experiments 
under  conditions  where  edge  diffusion  has  to  be  taken  into 
account.  The  influence  of  the  dimensionless  parameters*  $  and ®', 
the  initial  integration  stepwidth,  AT*0.  the  number  of 
collocation  pooints*  N*  and  the  kind  of  polynomial  used  to 
approximate  the  real  solution  are  investigated  and  discussed.  An 
additional  transformat  ion  of  the  radial  coordinate  is 
introduced.  The  numerical  values  for  a  large  number  of 
simulations  with  different  sets  of  the  above  parameters  are 
compared  to  reference  values. 
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INTRODUCTION 

In  part  1  1 L ]  of  this  series  we  derived  the  theory  for  the 
application  ot  orthogonal  collocation  to  the  s iauldtion  of  two- 
dimensional  models.  These  are  discribed  by  partial  differential 
equations  with  two  independent  space  variables  and  the  tine 
variable  t.  Let  us  consider  a  two-dimensional  diffusion 
situation  to  a  flat,  circular  electrode,  where  diffusion  not  only 
perpendicular  to  the  surface  may  occur.  This  additional 
diffusion  is  called  edge  diffusion.  The  concentration  c  of  each 
species  is  given  by  a  partial  differential  equation  of  the  form 

-5  *  D  fi-'  .  i  *£  .  “Jsl  -  Uci  111 
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where  s  is  the  distance  perpendicular  to  the  surface  and  r  is  the 
radial  distance  from  the  centre  of  the  electrode.  The  "kinetic 
term”  f(c>  is  determined  by  chemical  reactions  of  the  given 
species. 

Our  goal  is  to  calculate  c  as  a  function  of  x  and  r  at  any 
time  t  under  a  set  of  Initial  and  boundary  conditions.  The 
initial  conditions  model  the  experimental  starting  conditions, 
l.e.  in  general  only  the  starting  compound  is  present  in  a 
homogeneous  solution.  The  boundary  conditions  specify  the 
particulars  of  the  experiment  performed,  e.g.  chronoamperometry, 
cyclic  voltammetry  or  another  elect roans ly t ical  method. 

Orthogonal  collocation  12. )|  approximates  the  real  solution 
of  partial  differential  equations,  e.g.  ( 1 ) ,  by  orthogonal 
polynomials,  such  that  the  differential  equation  is  fulfilled  at 


certain  points,  the  collocation  points.  This  provides  a  method 

|  to  reduce  the  partial  differential  equations  to  e  system  ot 

I 

coupled  ordinary  differential  equations. 

In  part  1  ot  this  series  111  we  used  e  product  of  two 
polynomials  in  x  and  r  respectively  as  approximation  function  to 
the  real  solution  of  our  two •dimensional  model  for  edge 
effects.  We  solved  a  system  of  1  equations  of  type  (1) 
describing  an  EC  mechanism 

-ne  ~ 
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under  the  boundary  conditions  valid  Cor  chronoamperometric 
experiments. 

This  second  part  of  our  aeries  demonstrates  the  application 
of  the  derived  theory  to  chronoamperometry.  examine  the 
influence  of  several  parameters  on  the  numerical  solutions  of  the 
simulations!  the  dimensionless  parameters  «  ends' ,  the 
integration  stepwidth  *T'0«  the  number  of  collocation  points  N, 
and  the  type  of  polynomial  we  ere  using  to  approximate  the  real 
\  solution.  Also,  we  investigate  the  effect  ot  choosing  a 

|  different  dimensionless  coordinate  system. 

The  calculations  herein  have  been  performed  using  only  the 
concentration  of  the  substrate  A  at  the  collocation  points, 
bv  cause  the  current  through  the  electrode  may  be  calculated  from 
the  concentration  profiles  of  A  only.  The  concentration  profiles 
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o(  ipeclea  S  and  C,  Irowwt,  may  bn  calculated  frnm  the  [ornulae 
In  part  1  111 . 

He  coape  re  our  nuaerical  aiaulatlona  to  the  reaulta  derived 
by  an  aaact  aolution  of  a  related  problem  by  Oldham  |4|.  These 
nuabera  compare  very  veil  with  eolutiona  calculated  by  Aoki  and 
Oateryounq  1 51  and  Bmklh.n,  et  al.  [61,  all  using  different 
approachea  to  compute  valuea  for  the  current  at  a  finite  disc 
electrode  under  chronoamperoaetric  conditions  with  consideration 
of  edqe  effects.  All  these  approachea,  however,  are  restricted 
to  a  very  simple  mechanism,  the  reversible  electron  transfer. 

The  method  of  orthoqonal  collocation,  however,  is  -  once  derived 
for  a  specific  type  of  model,  ».q.  a  two-dimensional  model  -  very 
easily  used  to  simulate  all  possible  kinds  of  electrode 
mechanism.  Mt  only  have  to  change  the  kinetic  terns  f(cl  in  our 
Partial  differential  equations  II).  equations  describing  other 
species  than  A.  ■  or  c  are  very  similar  In  form  to  those  already 
derived  111.  This  will  be  shown  in  part  J  of  this  series  171  for 
the  boundary  conditions  of  cyclic  voltammetry.  Alao.  changes  in 
aechantm  require  only  minor  changes  of  the  computer  code  used 
to  calculate  numerical  data.  Thus,  orthogonal  collocation  proves 
to  be  a  flealble  method  td  generate  simulated  values  for 
electrochemical  eaperlaents  modelled  by  two-dimensional  diffusion 


The  Dirnews  lomLess  PAaMtereits  a  and  »• 

1t>  perform  the  method  of  orthogonal  collocation  we  have  to 
normalise  our  problem  Into  a  space  10. 1)  In  both  directions  s  and 


r.  Also,  we  have  to  divide  the  r*  space  into  an  inner  and  an 
outer  segment,  aie-iated  by  two  different  functions.  We  denote 
the  new  distance  coordinate  perpendicular  to  the  electrode  as  X« 
and  refer  to  the  normalized  inner  radial  se9ment  as  R  while  using 
R'  for  the  outer  coordinates.  The  dimensionless  simulation  time 


After  substituting  the  transformation  equations  given  in 
part  l  |1|  we  arrive  at  normalized  differential  equations  which 
contain  two  dimensionless  parameters,  •  and  t ' I 


Here^  rg  is  the  radius  of  the  electrode  disc*  imbedded  in  an 
insulator.  The  quantity  I*  is  a  distance  in  X  direction  from  the 
electrode  at  which  no  diffusion  occurs  during  the  simulated 
experiment*  while  M  is  an  analogous  distance  in  r-dlrectlon. 

Thus ,  L  and  M  define  a  region  in  real  space  which  is  simulated  in 
the  calculations.  The  reel  points  corresponding  to  our  N 
collocation  points  are  distributed  within  this  space.  Mr 
consider  two  limiting  cases  for  L.  An  analogous  discussion  may 
be  eppl  led  to  N. 

first*  L  may  be  very  large.  Then  the  points  corresponding 
to  the  collocation  points  in  real  apace*  l.e.  those  points  at 
which  the  concentration  of  the  species  involved  is  calculated, 
arc  spread  over  a  large  distance.  The  diffusion  layer  in  which 
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the  concent  ret  ion  of  a  species  is  depleted  ur  increased.  however, 
is  very  narrow  in  the  beginning  ot  the  simulation.  It  increases 
with  t*/^.  Thus,  at  short  times,  cor  respond  in- j  to  small  values 
of  the  dimensionless  simulation  time  T* ,  the  diffusion  may  have 
reached  only  few  or  none  ot  the  points  at  which  the  concentrat  ion 
is  c«.  ulated.  Consequently,  the  numerical  values  calculated  may 
not  represent  s  good  approximation  to  the  real  solution  of  our 
model  and  yield  erroneous  tsaults. 

On  the  other  hand,  if  L  is  very  large,  we  are  looking  at  a 
very  narrow  layer  of  solution  Ln  real  space.  Even  after  a  short 
time  the  diffusion  layer  may  have  reached  all  of  our  points  and 
the  boundary  condition  that  no  diffusion  occurs  at  x  >  1 
(corresponding  to  i  *  LI  is  no  longer  true.  We  are  using  the 
Mathematical  formulation  of  this  boundary  condition  (c£  (X-l)  • 

II  in  our  calculations  and  thus  the  numerical  values  may  deviate 
from  the  real  solution  after  this  is  happening. 

He  expect  that  the  approximation  solution  given  by 
orthogonal  collocation  is  valid  only  in  a  certain  time  range  and 
that  this  range  depends  on  the  values  of  t  and  $•  chosen. 

For  all  calculations  referred  to  in  ^his  paper  we  used  the 
condition  I  •  I*.  He  emulate  a  square  space  adjacent  to  the 
insulator,  because 

*  ,  <"  -  r0'2 

-  ■  1  ■  - —  (61 

or 
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L  “  ( H  -  rQ>  (7) 

Figures  la-c  show  calculated  values  for  the  current  obtained 
including  edge  effects  divided  by  the  ideal  Cottrell  current. 
Comparing  simulations  using  Legendre  polynomials  of  degree  N  •  6 
in  both  dimensions  but  different  values  for  1*0*  with  the 
reference  curve  given  by  Oldham  141,  we  see  that  the  above 
analysis  is  correct.  For  6*1  good  agreement  over  a  wide  range 
of  the  dimensionless  simulation  time  is  found  (2  ■  10"2  <  Ta  < 

1.1  «  10*l>.  To  give  a  specific  example  (D  •  1  *  10"*  cm2/s,  rQ 
*  0.1  cm I,  this  range  corresponds  to  real  times  t  between  2«10^ 
and  l.lxlO2  s  according  to  the  transformation  equation  relating  t 
and  T*  |1).  At  short  times  the  error  is  considerable  and  it  can 

be  seen  that  the  solution  needs  a  number  of  integration  steps  to 

become  stable.  At  *onger  times  the  collocation  simulation  raault 
increases  faster  than  Oldham’s  curve.  This  simulation  needed  40 
s  epu-time  to  calculate  the  entire  curve  (10**  <  T4  <  l.i"10*2l« 

More  accurate  simulations  for  shorter  times  may  be  obtained 
using  larger  values  for  For  example,  with  $  ■  10  in  a  range 

u:>tween  T’  •  2  ■  10*  2  and  1.2  »  10"2  good  agreement  with  Oldham’s 

curve  is  found,  while  |  ■  100  gives  a  reasonable  approximation 
for  2  >  10*4  4  t*  4  1.1  >  10"J.  In  the  latter  cases  significant 
deviation  from  Oldham’s  curve  is  seen  et  a  time  when  our 
numerical  values  increase.  This  occurs  at  the  same  time  as  the 
calculated  concentration  at  one  of  the  outmost  collocation  point 
in  the  x-direction  becomes  smaller  than  about  0.98,  thua 
confirming  our  conclusions  drown  above.  With  increasing  $, 
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however*  the  oncitlatinq  deviation*  at  the  beginning  occur  m  a 
narrower  time  frame.  again  in  accordance  with  our  analyst*. 

The  saw  behavior  nay  be  seen  it  we  use  polynomials  with  a 
larger  number  of  collocation  points  (N  •  V  12»15j  Figure*  2- 
4).  Now  the  "critical  concentration*  at  the  outmost  collocation 
points  moves  nearer  to  unity,  being  about  0.946  for  n  •  4.  The 
time  T'  at  which  the  deviations  occur,  however,  is  almost  the 
Sam  as  for  N  •  i.  This  leads  to  the  empirical  conclusion,  that 
the  upper  1 1 me  Halt  of  the  simulations  is  only  determined 
by  a.  An  analogous  approch  has  been  used  by  Whiting  and  Carr 
121  •  They  defined  the  maximum  simulation  time  which  is 
eccesslble  with  a  certain  •  as  T»  «  0.1 5/ a  .  In  our  case  we 
have  the  condition 

T’mas  *  *•*/•  <•> 

to  be  within  the  usable  time  range. 

THE  NUMBER  Of  COLLOCATION  POINTS  N 

In  the  preceding  paragraph  we  have  seen  that  the  number  of 
collocation  points  N  has  a  marked  influence  on  the  simulated 
numerical  values.  Thus,  we  compare  curves  obtained  with 
different  numbers  N  but  using  the  same  other  parameters  in  Figure 
9. 

As  N  increases,  the  net  of  collocation  points  becomes  more 
narrow,  forcing  the  approximation  function  to  fulfill  the 
differential  equation  at  more  and  store  points  exactly.  This  has 


effects:  first,  the  overall  solution  becomes  more  accurate 

but  w  also  have  more  simulation  pointy  near  the  boundaries 
fusing  Legendre  polynomials,  the  collocation  points  are  more 
densly  spaced  near  X  •  0  and  x  •  1  than  around  X  -  0.5).  Thus, 
for  short  times,  we  have  more  points  to  represent  the  steep 
concentration  profile  exhibited  shortly  after  applying  the 
potential  pulse  and  we  expect  the  range  within  which  the 
oscillating  deviations  from  the  reference  solution  occur  to  be 
more  narrow.  With  increasing  n,  the  overall  solutions  should 
reach  s  common  limiting  value  which  should  be  a  good 
approximation  to  the  real  solution. 

Comparing  the  results  given  In  Figure  5  we  see  that  this  is 
true.  By  increasing  the  number  of  collocation  points  from  N  ■  ( 
to  N  •  15  the  range  in  which  no  useful  results  are  to  be  obtained 
because  oscillations  occur  decreases  by  s  factor  of  about  10* 
giving  very  good  approximations  from  2  *  10“3  -  J  ’  10~2 
dimensionless  time  units.  The  upper  limit  of  the  useful  range  is 
not  affected,  as  Figures  1-4  show. 

Also,  the  results  differ  only  little  0.4%  when  increasing 
N  from  12  to  15)  with  increasing  N.  This  shows  the  consistency 
of  the  calculations.  The  "limiting"  values  reached  are  very 
close  to  Oldham's  solution  | 41 . 

Increasing  the  number  of  collocation  points  gives  more  and 
more  accurate  results  as  is  to  be  expected.  However,  the  epu 
time  used  to  perform  the  simulations  in  a  given  range  also 
increases  considerably.  While  integrating  a  system  o£ 

2  ■  6  ■  6  ■  72  ordinary  differential  equations  In  the  case  of  N  • 
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6  we  have  to  deal  with  162.  288  and  450  such  equations  it  we  use 
9.  12  or  IS  interior  collocation  points,  respectively.  The 
average  cpu  tiee  needed  to  calculate  one  concentration  protile 
increases  by  4  factor  of  10  when  N  is  increased  fro«i  6  to  15. 

THE  INTEGRATION  STEPWIDTH 

The  nest  parameter  we  want  to  discuss  is  the  initial 
integration  stepwidth,  AT'Q.  This  is  an  input  parameter  for  the 
subroutine  performing  the  integrations  of  our  systee  of  ordinary 
differential  equations.  The  value  choosen  is  used  by  the  , 

subroutine  to  calculate  the  first  concentration  profiles  I 

at  AT*^.  The  subroutine  DHpCGL  (in  the  SSP  library  package)  used 
by  us  then  “decides"  if  sufficient  eccurecy  is  available  with  | 

this  value.  If  not,  the  stepwidth  is  divided  by  2  and  the  J 

calculation  is  repeated.  Other  Integration  subroutines,  e.g.  | 

STIPF3  (3|,  also  have  the  option  to  increase  the  initial  j 

stepwidth  if  the  error  ie  below  e  specific  margin.  j 

To  evaluate  the  Influence  of  AT'Q  0(1  the  accuracy  and  ! 

performance  of  the  edge  effect  calculations  we  run  simulations  j 

with  N  •  i  and  9  using  it*0  values  ranging  from  10*?*  to  10*5.  ! 

Results  for  $  »  $•  •  1  ar«  given  in  Tables  1  and  2.  Me  see  | 

immediately  that  the  actual  value  of  AT'0  has  only  a  very  small  I 

effect  on  the  numerical  simulation  results.  Deviations  are  J 

0.0000?  in  the  worst  case,  this  corresponds  to  0.007t.  However, 
for  e  smell  iT'a  < 10“*  dlmenslonlese  time  unite),  the  cpu-time  ! 

needed  to  simulate  the  model  up  to  <  certain  time  increases  j 

considerably.  The  savings  using  AT*0  values  >  10" 3  are  only 


12 

about  10-201.  With  large  initial  integration  stepwidth*  there  is 

not  enough  accuracy  at  the  beginning  of  the  simulation,  thus 

causing  DHPCGL  to  cut  down  AT'0  to  smaller  Values.  The  actual 

choice  of  AT'  ,  however,  seems  not  to  be  critical, 
o 

COMPLETE  SIMULATIONS  USING  LEGENDRE  POLYNOMIALS 

After  discussing  this*  of  the  parameters  used  in  simulations 
of  chronoaaperometric  experiments  with  cons ldsrst ion  of  edge 
effects,  we  may  now  choose  sets  of  parameters  to  obtain  sn 
approximate  solution  of  the  differential  equations  over  a  large 
range  of  the  dimensionless  simulation  time.  The  number  of 
interior  collocation  points  will  be  15.  We  will  piece  together 
the  curve  from  solutions  calculated  with  different  ft's,  choosing 
the  appropriate  initial  integration  stepwidth.  The  trial 
functions  will  be  imgcndre  polynomials. 

The  result  is  given  in  Figure  6,  where  the  calculated 
current  divided  by  the  "ideal"  Cottrell  current  is  plotted 
against  T*  (in  logsrlthmic  seals).  Four  different  values 
for  K1000,  100,  10,  1)  have  been  used  and  the  entire  curve 
comperes  very  well  with  Oldham's  results  (4).  fti*  total  cpu-time 
needed  to  calculate  our  curve  was  about  1340  s  on  an  Amdahl 
470/V6  computer.  We  cover,  however,  more  then  4  orders  of 
magnitude  of  the  dimensionless  simulation  time.  Using  the  same 
constants  tor  D  and  r^  as  above,  we  simulated  the  response  of  sn 
electrode  with  edge  effects  from  4  ms  up  to  30  s.  At  the  lower 
limit  (T*  «  2.5»10~*)  the  current  with  edge  effects  deviates  only 
by  about  0.05%  from  the  Cottrell  current,  while  at  the  upper 
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Unit  (T*  «»  3*10“2)  the  current  in  real  experiments  may  be 
affected  by  convection,  thus  givinq  rise  to  anothe*-  kind  of  non¬ 
ideal  behavior.  Therefore,  we  can  caver  the  entire  accessible 
time  range  In  which  edge  effects  contribute  to  the  ex|>er imental 
current  with  this  simulation. 

OTHER  POLYNOMIALS 

Orthogonal  collocation  may  be  performed  using  a  variety  of 
trial  functions  as  long  as  we  have  a  set  of  complete  and  linearly  , 

independent  functions  (S| .  A  certain  class  of  sets  fulfilling  ' 

these  conditions  are  the  Jacobi  polynomials  p*5,c)  (x)  given  by  . 

Ml  I 

Pi4**’!*!  •  ?  (-I)**1***4.  y0  -  1  (9) 

"  1-0  i  0 

with  the  orthogonality  condition 

r)U)PNt*)  -  0  j  •  0,  1...,  M-l  (10) 

The  two  terms  »c  and  (1-x)4  are  “welght^actors*  which  determine 

the  behavior  of  the  function  at  the  boundaries.  For  c  •  0  the 

function  goes  to  xero  at  x  -  0.  for  4*0  sero  is  reached  for  x  »  . 

1.  In  our  case  the  radial  functions  have  nonzero  values  at  both  j 

boundaries.  Also,  the  x-functions  in  our  outer  space  have  a 

finite  value  at  a  ■  0  (the  insulator). 

Thus,  ms  would  expect  the  best  set  of  trial  functions  to  be 
Legendre  polynomials,  with  s  -  c  •  0.  He  used  Legendre 


polynomials  in  all  previous  calculations  in  this  paper,  but  tried 
also  several  other  polynomials  for  the;  radial  discretization 
function  (e.g.  6*0,  c  ■  2  and  4  •  1»  c  *  0.5  which  yielded  good 
results  for  problems  with  spherical  symmetry). 

The  calculations  showed  that  in  all  these  cases  the 
simulation  results  gave  poorer  agreement  with  Oldham's 
solution.  Me  conclude  therefore  that  the  use  of  Jacobi  type 
functions  other  than  Legendre  polynomials  cannot  b«  recommended 
in  our  case,  due  to  the  form  of  the  boundary  value  problem. 


ADDITIONAL  TRANSFORMATION  OF  THE  RADIAL  VARIABLE  R 


After  having  simulated  the  concentration  profile  In  our 
model  we  wish  to  calculate  the  current  through  the  electrode  as  a 
function  of  time.  In  order  to  compare  this  current  to  the 
•ideal"  Cottrell  current  we  derived  the  expression  (11)  HI* 


_ 

it*/*' 


247F  r'*  d«  un 

•  »X  X-0 

Cottrell 

where  we  have  to  Integrate  the  product  of  R  and  the  flux  to  the 
electrode  surface  over  the  entire  radius  of  the  circular 
electrode. 

This  integration  may  be  performed  by  a  quadrature  formula 
Mr*2 

/'  y(t)  w.  (r)dr  •  I  a.  Y  (l24 

•  *  1-1  ^  • 


where  the  Qj  depend  on  the  trial  function  used  |1|. 

The  form  of  the  integral  In  equation  (11)  would  suggest  a 
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it 


trial  (unction  with  6  «  0  and  t  »  1  •  We  have  seen  in  the  last 
paragraph,  however,  that  t)ie  only  suitable  trial  tunctions  tor 
our  problem  are  Legendre  polynomials  with  4  -*  c  *  0.  Thus,  it  is 
desirable  to  bring  (11)  into  a  tom  with  6  -  c  «  0  by  a  variable 
transformation. 

This  may  be  done  by  the  substitution 


Tlte  result  tor  species  A  is 


*CA 

3T* 


a2  .* 

>«•  *  .ft* 


and  the  boundary  condition  changes  to 


(  it) 


Cl-  x1 


Then 


(13) 


(*ii\  .  'W»<*A 

V* Alt**  1  '***  '  l>'*0 


(171 


do  •  a 

zOt1/J 


d«J  -  ‘  d Cp 
«#. 


If  we  substitute  these  equations  for  equations  (26)  and  (SO) 
(1^)  in  part  l  ill  and  repeat  the  derivation  of  the  discretized 

formulae  (9ft)  and  (104)  in  part  l  | 1 )  we  find  a  new,  very  similar 
system  of  ordinary  differential  equations. 

In  fact,  we  only  have  to  change  the  definition  of  the  matns 
elements  Ej#k  to 


The  integral  m  equation  <li)  reduces  to 


1 

2 


,i‘  i-c-) 

0  IX  X-0 


d  A 


thus  havinq  .  (or.  consistent  with  a  ■  c  -  0. 

H.  hava  to  apply  t ran. Corn. t ions  (U)-(1M  to  the 
differential  aquation  dascrlbinq  the  Inner  part  of  our  Kxiel  (see 
equation  (2t>.  part  1  of  this  aeries  111)  and  to  the  boundary 
condition  «t  I  ■  |  I X*  ■  0),  where  we  place  toqether  the  inner 
end  outer  function  (equation  (50)  in  part  1  111). 


i 

I 

I 

I 

I 

I 

I 


Bl.k  *  2-Cl,k  ^  4^Dl,k  ^ 

and  substitute  / ft V 2  for  /fF  in  the  definitions  of  the  vector 
elements  Gj,  the  Mtrix  elements  Li,*  and  the  discretized 
expression  (98). 

We  calculated  numerical  values  for  this  system.  The  results 
are  given  and  compared  to  the  nontrensformed  system  in  Figures  7- 
9  for  n  •  6,  9,  12.  Values  for  N  •  IS  are  compared  in  Itoble  3, 
while  Figure  10  gives  a  curve  pieced  together  using  4 
calculations  performed  with  N  *  15  and  different  ft  values.  The 


i 


I  - 

/ 


t? 
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numerical  values  are  consistent  with  changes  <  0.21  with  N 
increasing  by  3  In  a  region,  where  the  oscillatory  behavior  at 
the  beginning  has  stopped.  Coopering  to  the  nontransformed 
system  and  to  Oldham's  reference  values,  we  find  that  the  new 
equations  give  better  values  for  a  given  number  of  collocation 
points.  For  N  «  15  deviation  from  the  reference  values  is  in  the 
range  of  0.161  (Table  3). 

Figure  10  also  shows  how  well  our  numerical  values  compare 
to  Oldham's  model  over  more  than  4  orders  of  magnitude  in  T' . 

While  the  nontransformed  equations  gave  rise  to  a  curve  which  was 
lower  than  Oldham's  at  shorter  times  and  higher  at  longer  times 
(Figure  6),  we  see  that  our  transformed  system  gives  values  which 
are  slightly  lower  than  the  reference  over  the  entire  range. 

This  example  clearly  shows  that  the  choice  of  trial 
functions  and  transformations  has  an  influence  on  the  exact 
numerical  solutions  calculated.  Good  approximations,  however, 
may  be  available  with  the  "standard*  method  also  (i.e.  without 
consideration  of  the  type  of  functions  used).  A  more  detailed 
Investigation  on  the  effect  of  using  different  sets  of  trial 
functions  is  currently  under  way.  *** 

CONCLUSIONS 

Our  two~dimensional  diffusion  mode  1  using  orthogonal 
collocation  has  proved  to  provide  good  approximations  for  the 

I 

current  expected  at  a  flat,  circular  electrode  imbedded  into  an  j 

insulator  over  several  orders  of  magnitude  for  the  dimens  ionless  I 

time  variable  T* .  discussion  of  parameters  in  this  paper 

i 


shows  that  even  with  only  6  interior  collocation  points,  l.e.  a 
grid  of  2  ■  6  «  6  •  72  collocation  poUi\£  In  the  entire 
simulation  space,  reliable  nuaterical  values  may  be  obtained. 

With  increasing  number  of  collocation  points  the  simulations  are 
consistent  and  give  even  better  approximations.  However,  longer 
cpu-tlmes  and  costs  to  perform  the  calculations  are  necessary. 

We  again  like  to  point  out  that  changes  to  different  mechanisms 
requires  only  the  exchange  or  addition  of  a  few  statements  in  the 
computer  code.  Thus,  orthogonal  collocation  may  be  recommended 
as  a  fast  and  powerful  method  to  simulate  edge  effects  at 
electrodes. 


I  • 
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Table  1.  Numerical  values  for  simulation  of  a  chronoamperometric 
experiment  with  edge  effects:  current  with  edge 
effects  divided  by  Cottrell  current  simulation 
parameters:  N  »  8,  0  •  0*  •  1,  Lr-  -  polynomials  in 
both  dimensions,  AT*0  as  indicav*< 


TV«T-0 

10-2 

io-J 

10_s 

0.01 

1.1)249* 

1.132*6 

1.13256 

0.02 

1.2)991* 

1.23884 

1.23886 

0.03 

1.29917* 

1.29816* 

1 

1.29817 

0.04 

1.14796* 

1. 34787* 

1. 34787 

o.os 

1.19292* 

1. 39292 

1.39292 

~ 

0.06 

1.4)441* 

1.43444 

1.43444 

- 

0.07 

1.47)09 

1.47303* 

1.47305 

0.08 

1.90926 

1. 50926 

1.50926 

0.09 

1.94)49 

1 . 54 )S0* 

1.54349 

" 

0.1 

1.97606 

1.57607 

1.57607 

— 

*  Values 

calculated  by 

interpolation. 
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Tabic  2.  Numerical  values  for  simulation  of  a  chronoanperomet  r ic 
experiment  with  edge  effects:  current  with  edge 
effects  divided  by  Cottrell  current;  simulation 
parameters:  N  •  9,  &  *  $'  »1 «  Legendre  polynomials  in 
both  dimensions,  AT'0  as  indicated. 


TVfiT*0 

10"2 

10-3 

10-4 

10'5 

0.01 

1.17081* 

1.17080 

1. 17082 

1.17082 

0.02 

1.24847* 

1.24848* 

1.24848 

— 

0.0) 

1.30911 

1.30911* 

1.30911 

-- 

0.04 

1. 36075 

1.  36075* 

1. 36076 

— 

0.05 

1.40674 

1.40672 

1.40673 

- 

0.06 

1.44862* 

1.44864 

1.44864 

— 

0.07 

1.48748 

1.48748* 

1.48748 

- 

0.00 

1.52387 

1.52)87 

1.52)87 

— 

0.09 

1.55826* 

1.55826 

1.55826 

- 

0.1 

1.59100 

1.59100* 

1.59100 

- 

*  Va  1  ues 

calculated  by 

Interpolation. 

Table  1.  Numerical  values  for  simulation  of  a  chronoamperometr ic 
experiment  with  edge  effects;  current  with  edge  effects 
dev idi*d  by  Cottrell  current;  transformed  simulation  model  vs 
non transformed  system,  s  •  is ,*-»•-  10 ,  Legendre 
polynomials  in  both  dimensions. 


T* 

Oldham 'ssoiut ion  14) 

transformed  node  1 

non transformed  model 

»• 

10'4 

1.0175 

1.0206 

1.0192 

2 

1 .0249 

1.0221 

1.0210 

3 

1.0)05 

1.0291 

1.0270 

4 

1.0)61 

1.01)8 

1.0126 

5 

1.0)95 

1.0100 

1.0)67 

6 

1.04)1 

1.0416 

1.0405 

7 

1.0468 

1.0462 

1.0441 

0 

1.0500 

1.0406 

1.047) 

9 

1.05)1 

1.0616 

1.0504 

1. 

10-3 

1.0560 

1.0644 

1.05)4 

2 

1.0792 

1.0776 

1.0770 

3 

1.0970 

1.0964 

1.096) 

1  ' 
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fi cuke  Captions 

Figure  1.  SmuUUona  ot  chronoamperometr ic  e»|>enaenta 

including  edge  tttecls  with  orthogonal  collocation; 
smalation  parameters;  number  ut  collocation  joints 
in  each  diMnsion  N  *  6,  initial  integration 
stepwidth  at'0  •  1.10~4,  Legendre  polynomials  m  both 

dimensions,  ( - )  a  »  100,  ( - )  »  >  10.  < - )  c 

*  1;  (....)  values  calculated  according  to  Oldham 
l4l ■  The  arrows  indicate  T* ,  when  the  concentration 
at  one  ot  the  outmost  collocation  points  reaches 
0.98. 

Figure  2.  Simulations  ot  chronoampe rone t c ic  exper i merits 

including  edge  effects  with  orthogonal  collocation; 
simulation  parameters:  number  of  collocation  points 
in  each  dimension  N  »  9,  other  parameters  as  in 
Pigure  1.  The  arrow  indicates  T' ,  when  the 
concentration  at  one  ot  the  outmost  collocation 
points  reaches  0.996. 

Figure  3.  Simulations  of  Chronoampe rome trie  experiments 

including  edge  effects  with  orthogonal  collocation; 
simulation  parameters;  number  of  collocation  points 
in  each  dimension  N  »  12 ,  other  parameters  as  in 
Figure  1. 

Figure  4.  Simulations  of  chronoampe rome tr ic  experiments 

including  edge  effects  with  orthogonal  collocation; 
simulation  parameters;  number  ot  collocation  points 
in  each  dimension  n  »  IS,  other  parameters  as  in 


24 

Figure  1.  The  ordinate  scale  is  expanded  compared  to 
tl>e  other  figures. 

Figure  5.  Simulations  of  chronoampe rome t r ic  experiments 

including  edge  effects  with  orthogonal  collocation; 

simulation  parameters;  ( - )  M  -  6,  (—'—■•—)  N*9.  t  - 

- — )  M  -  12,  N  ■  IS;  initial  integration 

stepwidth  AT ‘ Q  •  1 . 10“4  <  Legendre  polynomials  in  both 
dimensions,  dimensionless  parameters  $  -  e*  -  1, 
(....)  values  calculated  according  to  Oldham  (4|. 

Figure  6.  Simulations  to  cover  the  entire  accessible  range  ot  t 

tor  an  example  with  D  -  10"*  cm2/s  and  rg  •  0>1  cm. 

pieced  together  from  4  calculations  with  N  •  IS  and 

different  0  compared  to  Oldham's  solution  (••*); 

points  calculated  with  0  •  »•  >1000  (••*),  0  ■  0*  * 
/!«** 

1MIOOO)  -10  (QQQ)  ,  0  •  0*  -  1  (AAAI- 

Figure  7.  Comparison  of  simulations  using  the  model  system 

given  in  part  1  ot  this  series  111  and  the  formulae 
derived  from  the  transformation  <131/  H  * 

6,  0  -  0'  »  10;  (— — )  before  transformation,  t'- - ) 

after  transformat  ion ,  (....)  Oldham's  curve. 

Figure  8.  Comparison  of  simulations  using  the  model  system 

given  in  part  1  of  this  series  | 1 i  and  the  formulae 
derived  from  the  transformation  (13);  N  - 
9,  0  •  0*  >10;  ( — •—  )  before  transformat  ion, 
after  transf or mat  ion,  (,...)  Oldham's  curve. 
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CaabnARO,  Haaaachuaotta  0213* 
Larry  I.  Plot* 

Pavol  Uoaaoas  Support  Coat or 
Coda  30734.  Buildup  2*04 
Crono,  Indiana  47322 

S.  ftubv 

DOC  (STM) 

400  C  Stroat 
Vaahio*toa,  D.C.  20343 

Dr.  Aaroa  UaU 

Brown  Univaraity 

Dopartaaat  af  Chaaiatry 

Pro*  idonca,  Skoda  Inland  021*2 

Dr.  I.  C.  Chudacoh 
Me  Craw  Crfi  non  Camp  t Of 
fdiaon  Battary  Diviaioa 
Poat  Offica  ftoa  20 
ftlaaafiatd.  Paw  Jaraay  07003 

Dr.  A.  3.  lard 
Uaivoraif y  of  Tomoo 
Dopartaaat  af  Chaaiatry 
Auatia,  Taaaa  70712 

Dr.  H.  P.  Picholaaa 
flactroaica  Roaaarck  Caatar 
Bock  wall  latarnat ieael 
337*  Niraloaa  Avaauo 
Aaafaaia,  Califaraia 

Dr.  Dana  I'd  P-  Em»t 

Naval  Surfaca  Uaopono  Coat or 

Coda  ft- 33 

Uhi to  Oak  Laboratory 

Silvor  Sarin*,  Nary load  20*10 


Dr.  ft.  P.  Van  Dayan 
Do  part no at  of  ChOBkitrv 
Northaaotara  Dnivaraity 
Bvaaataa,  Illinois  40201 

h.  I.  ftaalnJftu 
DtMhM^T  Chaaiatry 
Uni**yrUyof~~ALh*£l« 

Alborta 
CANADA  fftC  2C2 

Dr.  Picbaol  J.  Paavar 
Dopartaaat  of  Chaaiatry 
Nichigao  Stata  Paivaraity 
Coot  Uaaia*.  Michipa*  40*24 

Or.  B.  David  ftauh 

1XC  Corporation 

33  Chaaol  Stroat 

Nawton,  Naaoachaaatco  02130 

Dr.  J.  David  Narpana 
Baoaarcb  Labor at or lao  Diviaioa 
Nucha*  Aircraft  coapaay 
1011  Malibu  Caayoa  toad 
Malibu,  Califaraia  *0343 

Dr.  Marria  Pltitrtota* 
Dopartaaat  af  Chaaiatry 
Paivaraity  af  SduthMptM 
hnuapita  30ft  W  b|lMd 

Dr.  Jaaot  Oat ary ana* 

Da part aa at  af  Chaaiatry 
Stata  Paivaraity  of  Pm 
fork  at  Buffala 
Buffalo,  Pau  York  14214 

Dr.  I.  4.  Oatavyoaa* 
Dopartaaat  af  Chaaiatry 
Stata  Paivaraity  af  (•» 

Tort  at  toff ala 
•uffala.  Pan  Tort  14214 

nr.  Jaaaa  ft.  NoOaa 
Naval  Uadamatat  fyataaa 
Caatar 
Coda  3432 

Nawpert ,  Ihada  la land  02*40 


j 
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TtcmuM.  MfonT  Bimmmox  list,  nt 


4?l:CAH:7|4:l.b 

7(1*4  77-909 


Nr.  I.  Nowak 

Naval  Jtaaaarc*  Laboratory 
Coda  *130 

Woaklagtoa.  D-C.  2017) 

Nr.  John  7.  RoallhM 
Skanaago  Tallay  Ca«*aa 
Faaaaylvanta  Stata  Dnlvaralty 
9 karoo,  fmnirltanU  1*14* 

Nr.  D.  r.  SkrSaar 
fttpanaat  of  Ckariltry 
Nortkwascara  Dnlvaralty 
Ifamtia,  Tlllaola  *0201 

Nr.  D.  N.  ukltaara 
Oarartaaat  of  Notarial*  klaaca 
Nortkwaatarn  Dnlvaralty 
Cvaoataa,  llllaoia  *0201 

Nr.  Alan  tavlck 
Naaartaant  of  Cka^atry 
Tka  Dnlvaralty 
Sontkaaptoa,  SO*  MR  England 

Dr.  A.  Stay 
•AVSCA-543) 

NC  #4 

2341  Jaffaraaa  Darla  Ntgkway 
Ar llAftmi,  Vlrglola  203*2 

Nr.  laka  KlncatA 
Nayartaaat  af  tka  Navy 
Stataglc  tyataaa  Trnjact  Offlea 
>008  *01 

We-klngton,  DC  2037* 

N.  L.  Kotor toon 
Nano gar ,  five  tracking  cal 
Naaar  Sonlcaa  Dtvlotaa 
Naval  Waayana  Soy port  Cantor 
Crana,  Indiana  47321 

Nr.  titan  Catrno 
tnargy  *  Envirannant  MoUlaa 
Uar**et  Bark  a  lay  Laboratory 
Dnlvaralty  af  California 
Bart* lay,  California  *4720 


Satis* 


Nr.  Barnard  Splalvogal 
U.S.  A ray  ftaaoarck  Offlea 
7.0.  Boa  122 II 

1  Raaoarck  TrUagl#  7ark,  NC  27709  1 

l*r.  Canton  Elliott 
Air  Forca  Offlea  of 
S<i*neltte  Raaaarck 
1  Rolling  ATI 

Wa*Mng con,  NC  20332  1 


Dr.  Davfd  Altana 
CkonAatry  Dapartaaat 
(  Ranaaalaar  Felytacknlc  Tnatltato 

Troy,  NT  121*1  1 

Nr.  A.  7.  B.  Lavar 
Ckaadatry  D*part«»nt 
1  Tork  Unlvaralty 

Oovnavlao,  Ontario  M3JI73  1 


Hr.  Maori ca  7.  Horpky 
1  Naval  Kao  Syatana  Co— a* 

*3*32 

2221  Jaffaraaa  Davla  Rlgkway 
Arlington,  7A  203*0  1 

Dr.  Stool a lav  Styak 
1  Naval  Oeaan  Syataaa  Cantor 

Coda  *343 

San  Dlasa,  CA  *3132  1 

Dr.  Granary  Farrington 
Oayartnant  of  Hacartala  Selanca  * 

1  eogtaaarlag 

Dnlvaralty  af  Nnnaylvanta 
TkiladnlyM*,  7A  1*104  1 


Nr.  Bmco  Nana 

Dayartaant  af  CagtaaorlaR  A 
1  Appllad  Sclanra 

Dnlvaralty  af  Caltfamla 

too  Aagalaa,  CA  *0024  1 
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TtCHHlCAL  RETORT  PIST11BDTI0N  LIST.  33* 


Dr.  Rieka  ToaAiwlc* 
Depart  boi  of  Tkyalea 
Brooklyn  Collage 
Brooklyn,  NT  11210 

kr.  Laaoar  Blua 
Depart  But  af  Fkyatee 
Unlvaralty  of  Toarta  Rico 
Rto  Ftadraa.  7R  00*31 


1 


kr.  Joaepk  Cor Aon  U 
IBM  Coryoratlaa 
R33/2R1 

3 400  Cattle  Road 

San  Joaa,  CA  *31*3  1 

Nr.  » chart  Soanana 
Jot  kropnltfen  Laboratory 
California  laatltota  af  Taeknolagy 
Faaadana.  CA  *1103  1 
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